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ABSTRACT 


Climatological  temperature  and  salinity  fields  are  used  to  calculate  the  salinity 
contribution  to  density  and  dynamic  height  fields  in  the  Leeuwin  Current  System  (LCS). 
While  the  temperature  gradient  is  primarily  linear,  with  warmest  water  to  the  north,  the 
salinity  fields  are  spatially  inhomogenous.  A  comparison  of  density  fields,  calculated 
with  constant  and  variable  salinity,  shows  that,  off  Western  Australia,  the  density  field  is 
primarily  determined  by  temperature.  Off  Southern  Australia,  the  density  field  is 
dependent  on  warm  and  salty  (subtropical)  and  fresh  and  cold  (sub-Antarctic)  water 
masses.  While  the  dynamic  height  fields,  calculated  with  constant  and  variable  salinity, 
show  similar  flow  patterns  off  Western  Australia,  different  flow  patterns  are  found  off 
Southern  Australia. 

In  addition  to  the  analysis  of  climatological  fields,  a  primitive  equation  ocean 
model  is  used  to  investigate  the  role  of  salinity  in  the  formation  of  currents  and  eddies  in 
the  LCS.  Two  identical  ocean  models,  one  with  a  climatological  salinity  field  and  the 
other  with  no  horizontal  salinity  gradients,  are  run  and  compared  with  each  other.  Despite 
the  model  runs  being  initialized  with  similar  temperature  distributions,  there  are  relatively 
large  temperature  and  density  differences  in  the  Southern  Australian  region,  due  to  the 
advection  of  water  masses  by  the  Leeuwin  Current. 

Based  on  the  climatological  analyses  and  the  results  of  the  model  experiments,  it 
is  concluded  that,  descriptively  and  dynamically,  both  temperature  and  salinity  are 
essential  to  accurately  characterize  the  large-scale  circulation  of  the  LCS. 
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I.  INTRODUCTION 


At  a  given  pressure,  the  density  of  sea  water  depends  on  temperature  and  salinity. 
Except  for  high-latitude  regions  and  near  sources  and  sinks  of  salt,  temperature  is  usually 
considered  to  vary  more  than  salinity  in  the  upper  ocean  and  is  generally  considered  to  be 
the  major  source  of  density  variations  (Pickard  and  Emery,  1990). 

Recent  dynamic  height  analyses  by  Batteen  et  al.  (1995)  have  shown  that  the 
distribution  of  salinity  in  the  California  Current  System  (CCS)  can  be  important  in 
defining  the  large-scale  circulation  of  the  CCS.  In  particular,  the  mean  variability  of 
salinity  has  been  shown  to  be  responsible  for  a  significant  equatorward  component  along 
the  coast  of  California  and  a  strong  offshore  component  adjacent  to  Baja.  Batteen  et  al. 
(1995)  concluded  that,  descriptively  and  dynamically,  both  temperature  and  salinity  are 
essential  to  accurately  characterize  the  large-scale  structure  of  the  CCS. 

Besides  the  CCS,  other  eastern  boundary  current  (EBC)  regions,  such  as  the 
Leeuwin  Current  System  (LCS),  could  show  significant  changes  if  salinity  as  well  as 
temperature  effects  are  considered.  The  Leeuwin  Current  is  an  anomalous  EBC  which 
flows  poleward  (against  the  prevailing  wind  direction)  along  the  Western  Australian 
coast,  down  to  Cape  Leeuwin  (See  Figure  1  for  geographical  locations),  and  swings 
eastward  extending  as  far  east  as  the  Great  Australian  Bight  (e.g.,  Cresswell  and  Golding, 
1980).  There  is  general  agreement  (e.g.,  Godfrey  and  Ridgway,  1985)  that  the  Leeuwin 
Current  is  generated  by  a  meridional  pressure  gradient,  which  overwhelms  the  opposing 
wind  stress.  The  source  for  the  Leeuwin  Current  water  is  predominantly  geostrophic 
inflow  from  the  west  (e.g.,  McCreary  et  al,  1986;  Thompson,  1987)  and  is  augmented  by 
a  source  from  the  North  West  Shelf  (e.g.,  Gentilli,  1972),  possibly  having  its  origin  in  the 
Pacific  Ocean  (Godfrey  and  Ridgway,  1985). 

Since  the  Leeuwin  Current  is  a  density  driven  current,  both  temperature  and 
salinity  could  contribute  to  its  development.  While  the  temperature  gradient  is  primarily 
linear,  with  warmest  water  to  the  north,  the  salinity  distribution  is  spatially 
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inhomogeneous  (Figure  2).  Off  of  Western  Australia  low  salinity  water  is  present  to  the 
north  and  northwest,  high  salinity  to  the  southwest,  and  low  salinity  to  the  south. 
Depending  on  the  time  of  year,  these  water  masses  could  contribute  to  the  development 
of  the  Leeuwin  Current,  subsurface  currents,  and  eddies. 

Here  we  ask  the  question  whether  salinity  variability  is  an  important  consideration 
for  defining  the  large-scale  circulation  of  the  LCS.  To  quantify  the  importance  of  salinity 
in  the  LCS,  we  calculate  its  contribution  to  density  and  to  dynamic  height  in  the  LCS 
(Chapter  II).  After  incorporating  salinity  into  an  ocean  model  (Chapter  III,  Section  A), 
two  identical  ocean  models,  one  with  a  climatological  salinity  field  and  the  other  with  no 
temporal  or  spatial  salinity  gradients,  are  run  and  compared  with  each  other  (Chapter  III 
Section  B).  Although  the  results  with  variable  salinity  included  show  similar  features  to 
the  results  with  constant  salinity  off  Western  Australia,  there  are  important  quantitative 
differences  in  the  Cape  Leeuwin  and  Southern  Australian  regions.  The  low  and  high 
salinity  tags  in  the  LCS  are  found  to  be  particularly  useful  for  identifying  the 
contributions  of  the  different  water  masses  to  the  development  of  the  Leeuwin  Current 
and  eddies. 
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II.  ANALYSIS  OF  HYDROGRAPHIC  AND  DYNAMIC  HEIGHT  FIELDS 

A.  ANALYSIS  OF  ANNUAL  HYDROGRAPHIC  AND  DYNAMIC  HEIGHT 
FIELDS 

Two  major  water  masses  occur  for  all  or  part  of  the  year  off  Western  Australia 
(Figure  2):  a  warm  (tropical)  low  salinity  water  mass  from  the  North  West  Shelf  (e.g., 
Weaver  and  Middleton,  1989)  possibly  having  its  origins  in  the  Pacific  Ocean  (Godfrey  and 
Ridgway,  1985),  and  a  cooler  (subtropical),  higher  saline  water  mass  called  the  'Western 
Australian  Current'  which  is  influenced  by  the  presence  of  a  large  semi-permanent  trough 
extending  in  a  northeastward  direction  (Andrews,  1977). 

The  seas  just  south  of  Cape  Leeuwin  and  eastward  to  the  entrance  (near  Esperance 
in  Figure  2)  to  the  Great  Australian  Bight  are  influenced  by  several  water  types:  1) 
subtropical  water  from  west  of  Western  Australia  carried  in  a  branch  of  the  Western 
Australian  Current  in  summer  (Andrews,  1977),  2)  tropical  water  from  north  of  Western 
Australia  brought  by  the  Leeuwin  Current  in  autumn  and  winter  (Rochford,  1969;  Cresswell 
and  Golding,  1 980),  and  3)  local  Sub- Antarctic  water  characterized  by  its  coldness  and  low 
salinity  (Rochford,  1986). 

To  explore  the  water  masses  of  this  area,  the  climatological  data  sets  of  Levitus  et 
al.  (1994)  and  Levitus  and  Boyer  (1994)  were  employed.  The  data  are  gridded  with  one- 
degree  horizontal  resolution  at  standard  depth  levels.  These  data  sets,  which  include 
annual  and  monthly  mean  temperature  and  salinity  fields,  were  used  to  compute  annual 
and  seasonal  density  and  dynamic  height  fields. 

Before  presenting  the  annual  and  monthly  climatological  temperature,  salinity, 
density,  and  dynamic  height  fields,  we  describe  how  the  calculations  for  density  and 
dynamic  height  were  obtained.  Following  Batteen  et  al  (1995),  the  density  is  computed 
using  the  following  approximate  equation  of  state  for  seawater: 

p  =  p,[\-a{T-T,)  +  f3{S-S,)], 
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where  a  and  P  are  the  expansion  and  contraction  coefficients  for  temperature  T  and 
salinity  S,  respectively,  p  is  density,  Sq  is  the  reference  salinity,  and  Tq  is  the  reference 
temperature.  This  equation  can  estimate  the  effect  of  salinity  for  a  given  region  if  typical 
values  of  a,  P,  Tq,  and  Sq  are  used.  For  the  LCS  region,  a  is  chosen  to  be  2.5x10'^  (°K)'', 
P  is  7.5x10"^,  To  is  278.2  °K  and  Sq  is  35.5  psu  (practical  salinity  units). 

To  assess  the  effect  of  salinity  on  density  in  the  LCS,  we  calculate  density  in  the 
above  equation  for  variable  and  for  constant  salinity.  Calculating  density  for  variable 
salinity  means  that  the  density  in  the  above  equation  is  computed  at  a  given  level  with 
climatological  values  of  T  and  S.  Calculating  density  for  constant  salinity  means  that  the 
density  in  the  above  equation  has  S  replaced  by  Sq,  so  that  the  above  equation  reduces  to 

p  =  p,[\-a{T-T,)]. 

To  assess  the  effect  of  salinity  gradients  on  the  large-scale  ocean  circulation  of  the 
LCS,  dynamic  height  fields  are  calculated  with  and  without  the  inclusion  of  salinity 
variations.  The  dynamic  height  anomaly  AO  is  calculated  by  integrating  the  specific 
volume  anomaly  5  over  depth  and  comparing  it  with  AO35  syp  (where  the  subscripts  35.5, 
T,  and  P  stand  for  salinity  35.5  psu,  temperature,  and  pressure,  respectively),  calculated 
by  integrating  by.  Comparisons  between  horizontal  gradients  of  AO  and  AO35  syp  are  used 
to  provide  insight  concerning  the  importance  of  salinity  variations  to  the  flow  field.  The 
reference  level  used  in  this  study  is  1000  m  depth. 

The  climatological  annual  mean  temperature  fields  are  shown  at  10,  100,  and  200 
m  depth  in  Plate  la,  and  at  400  m  depth  in  Figure  3a.  At  all  depths  shown,  the 
temperature  distribution  is  primarily  zonal,  with  the  isotherms  running  roughly  east-west. 
Maximum  temperatures  are  in  the  north  and  decrease  to  the  south.  The  temperature 
differences  from  north  to  south  are  ~  1 1°C,  9°C,  7°C,  and  2°C  at  10,  100,  200,  and  400  m 
depth,  respectively. 

The  climatological  annual  mean  salinity  fields  are  shown  at  10,  100,  and  200  m 
depth  in  Plate  lb,  and  at  400  m  depth  in  Figure  3b.  The  common  features  of  the  10  and 
100  m  depth  plots  can  be  clearly  seen:  low  salinity  (tropical)  water  to  the  north  and 
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northwest,  high  salinity  (subtropical)  water  to  the  southwest,  and  low  salinity  (sub- 
Antarctic)  water  to  the  south.  The  maximum  salinity  for  10  m  depth  is  centered  at  ~  32°S, 
while  the  maximum  salinity  for  100  m  depth  is  centered  at  ~  28‘’S.  The  salinity  gradient 
weakens  from  10  to  100  m  depth.  The  35.5  salinity  contour  extends  all  along  Southern 
Australia  at  10  m  depth,  but  only  as  far  east  as  Albany  at  100  m  depth.  At  200  m  depth, 
the  salinity  pattern  is  similar  to  the  zonal  temperature  pattern  with  maximum  salinity  in 
the  north  and  minimum  salinity  in  the  south.  At  this  depth,  the  35.5  salinity  contour 
extends  only  as  far  south  as  Fremantle  (~  32°S).  At  400  m  depth,  the  salinity  difference  is 
only  0.4  psu  from  north  to  south. 

The  annual  density  fields  calculated  with  constant  (variable)  salinity  at  10,  100, 
and  200  m  depth  are  shown  in  Plate  Ic  (Plate  Id).  The  annual  density  fields  calculated 
with  constant  (variable)  salinity  at  400  m  depth  are  shown  in  Figure  3c  (Figure  3d).  A 
comparison  of  the  upper  levels  (i.e.,  10  and  100  m  depth)  in  Plate  Ic  and  Plate  Id  shows 
the  following.  Although,  off  Western  Australia  the  density  is  primarily  determined  by  the 
temperature  field,  the  density  fields  in  Plate  Id  are  generally  less  dense  than  those  in  Plate 
Ic  due  to  the  presence  of  the  low  salinity  (tropical)  water  mass.  South  of  ~  34°S,  the 
effect  of  salinity  on  density  becomes  significant.  In  particular,  off  Southern  Australia,  the 
density  is  influenced  by  both  the  warm  and  salty  (subtropical)  water  mass  and  the  fresh 
and  cold  (sub- Antarctic)  water  mass.  At  200  m  depth,  the  density  in  Plate  Id  is  generally 
greater  than  that  in  Plate  1  c  due  to  the  presence  of  high  salinity  throughout  the  region.  At 
400  m  depth,  the  density  field  in  both  figures  is  relatively  uniform  throughout  the  region. 

Figures  4  and  5  show  the  annual  dynamic  height  patterns  at  10,  100, 200,  and  400 
m  depth  for  constant  salinity,  i.e.,  AO355JP,  and  for  variable  salinity,  i.e.,  AO, 
respectively.  In  both  figures,  the  upper  level  (i.e.,  10  and  100  m  depth)  dynamic  height 
contours  tend  to  be  perpendicular  to  the  west  coast  of  Australia,  which  is  consistent  with 
onshore  geostrophic  flow  toward  Western  Australia.  Between  34°S  and  36°S,  while  the 
upper  level  dynamic  height  fields  in  Figure  4  tend  to  be  zonally  distributed,  those  in 
Figure  5  tend  to  vary  in  both  the  meridional  and  zonal  directions,  due  to  the  influence  of 
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the  salty  (subtropical)  water  mass  offshore.  The  lower  level  (i.e.,  200  and  400  m  depth) 
dynamic  height  patterns  in  both  figures  show  equatorward  flow  off  Western  Australia, 
which  is  consistent  with  the  depths  at  which  the  equatorward  imdercurrent  has  been 
observed  (e.g.,  Thompson,  1984). 

B.  ANALYSIS  OF  SEASONAL  HYDROGRAPHIC  AND  DYNAMIC  HEIGHT 
FIELDS 

Since  the  LCS  has  a  strong  seasonal  (e.g.,  the  Leeuwin  Current  is  strongest  from 
March  to  August)  as  well  as  annual  cycle,  the  monthly  mean  temperature  and  salinity 
climatological  fields  are  examined  to  evaluate  the  effects  of  the  seasonal  evolution  of 
salinity  on  the  density  and  dynamic  height  fields.  The  annual  mean  temperature  and 
salinity  fields  showed  that  the  water  mass  characteristics  of  the  LCS  are  very  similar  at  10 
and  100  m  depth.  Since  the  seasonal  cycle  has  essentially  the  same  vertical  pattern  in  the 
upper  layers,  we  present  only  the  10  m  depth  seasonal  fields  here. 

The  same  three  water  masses  described  earlier  are  discussed  here  to  help  identify 
how  salinity  seasonally  influences  the  density  and  dynamic  height  fields.  These  water 
masses  are:  the  warm  and  fresh  tropical  water  mass  off  northwestern  Australia,  the  warm 
and  salty  subtropical  water  mass  off  southwestern  Australia,  and  the  cold  and  fresh  water 
mass  off  Southern  Australia. 

Plate  2a  shows  the  monthly  mean  temperature  fields  in  February,  May,  August, 
and  November.  As  the  November  and  February  plots  show,  the  temperature  gradient  off 
Western  Australia  during  the  austral  spring  and  summer  has  a  45  degree  slope  from 
northeast  to  southwest,  which  provides  the  source  for  an  onshore  geostrophic  current  (i.e., 
the  Leeuwin  Current).  As  the  May  and  August  plots  show,  the  slope  becomes  even 
steeper  during  the  austral  fall  and  winter,  which  is  consistent  with  the  Leeuwin  Current 
being  strongest  during  this  period.  Off  Southern  Australia,  except  for  the  advection  of  the 
Leeuwin  Current  in  the  vicinity  of  Cape  Leeuwin  during  the  austral  fall  and  winter,  the 
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temperature  pattern  follows  a  typical  seasonal  cycle  of  warmest  (coldest)  temperature  off 
the  coast  in  summer  (winter). 

Plate  2b  shows  the  monthly  mean  salinity  fields  in  February,  May,  August,  and 
November.  In  February,  the  salinity  field  shows  a  zonal  distribution  of  fresh  water  in  the 
north,  which  becomes  saltier  poleward  to  ~  32°S.  There  is  an  offshore  region  of 
maximum  salinity  of  -  35.8  from  ~  28°S  to  35°S,  centered  at  -  32°S.  South  of  -  35'’S, 
fresher  water  is  found.  In  May,  the  salinity  gradient  has  a  45  degree  slope  off  Western 
Australia.  The  maximum  salinity  has  shifted  southward  to  -  34°S.  In  August,  the  fresh 
water  off  northwest  Australia  has  moved  farther  southward,  while  the  maximum  salinity 
has  moved  farther  offshore.  In  November,  off  northwest  Australia,  the  fresh  water  has 
been  replaced  by  more  saline  water.  The  maximum  salinity  has  retreated  equatorwaxd  to 
-31°S. 

The  density  fields  shown  in  Figures  6a  and  6b  were  calculated  with  constant  and 
variable  salinity,  respectively.  In  general,  throughout  the  year,  the  density  fields  off 
Western  and  Southern  Australia  in  Figure  6b  are  less  dense  than  those  in  Figure  6a,  due 
to  the  presence  of  low  salinity  (tropical  and  sub-Antarctic)  water  masses.  In  the 
southwestern  Australia  region,  the  density  fields  in  Figure  6b  are  more  dense  than  those 
in  Figure  6a  because  of  the  presence  of  the  salty  (sub-tropical)  water  mass. 

Figures  7a  and  7b  show  the  seasonal  dynamic  height  patterns  for  constant  and 
variable  salinity,  respectively.  In  both  figures,  throughout  the  year,  off  Western  Australia, 
the  dynamic  height  flow  pattern  is  consistent  with  onshore  geostrophic  flow  toward 
Western  Australia.  In  contrast,  south  of  36°S,  while  the  dynamic  height  fields  in  Figure 
7a  tend  to  be  zonally  distributed  throughout  the  year,  those  in  Figure  7b  tend  to  be  vary  in 
both  the  meridional  and  zonal  directions,  due  to  the  influence  of  the  salty  (subtropical) 
water  mass  offshore. 

To  quantify  the  seasonal  effect  of  salinity  on  the  dynamic  height  field,  the 
dynamic  height  field  at  several  offshore  and  neairshore  locations  were  calculated.  Figure  8 
shows  the  meridionally  averaged  dynamic  height  field  at  108.5°E  from  1)  22°S  to  29% 
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2)  30°S  to  35°S,  and  3)  36°S  to  40°S.  In  region  1,  the  dynamic  height  with  variable 
salinity  (solid  line)  is  similar  to  that  with  constant  salinity  (dashed  line).  In  the  other  two 
regions,  the  dynamic  height  fields  are  dependent  on  both  temperature  and  salinity  from  ~ 
May  to  August  in  region  2  and  for  most  of  the  year  in  region  3.  In  region  2,  the  seasonal 
dependence  on  both  temperature  and  salinity  coincides  with  the  time  frame  that  the 
Leeuwin  Current  is  strongest.  In  region  3,  the  seasonal  salinity  effect  is  consistent  with 
the  time  that  the  subtropical  salty  water  moves  poleward  to  ~  40°S. 

The  dynamic  height  fields  calculated  at  several  nearshore  locations  are  shown  in 
Figure  9.  In  the  nearshore  regions  1  and  2,  which  correspond  to  coastal  locations  off 
Western  Australia,  because  the  salinity  gradient  is  not  as  strong  as  in  the  offshore  regions, 
the  salinity  effect  is  not  very  significant  in  these  regions.  In  contrast,  in  the  nearshore 
region  3,  which  corresponds  to  a  coastal  location  off  Southern  Australia,  the  dynamic 
height  is  dependent  on  both  temperature  and  salinity  throughout  most  of  the  year, 
indicating  that  in  this  region  salinity  plays  a  very  important  role. 
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III.  EFFECT  OF  SALINITY  ON  AN  OCEAN  MODEL  OF  THE  LCS 

A.  THE  MODEL 

1.  Model  Description 

The  numerical  model  in  this  study  was  originally  used  for  a  coarse  resolution, 
closed  basin  by  Haney  (1974),  and  later  adapted  for  eddy-resolving,  limited  EBC  regions 
with  open  borders  on  the  northern,  western,  eastern,  and  southern  boundaries  by  Batteen 
et  al.  (1989,  1996)  and  Batteen  (1996).  The  limited  area  EBC  model  is  multi-level,  uses 
non-adiabatic  primitive  equations  on  a  beta-plane,  and  has  both  baroclinic  and  barotropic 
velocity  components.  The  model  is  based  on  the  hydrostatic,  Boussinesq,  and  rigid  lid 
The  governing  equations  are  as  follows: 


du 

dt 

(1) 

dv 

dt 

(2) 

o 

ii 

(3) 

1 

11 

(4) 

P  =  Po{\-a{T -T,)+ P{S -S,)] 

(5) 

dT  , 

dt  ^  ^  dz^ 

(6) 

dS  , 

—  =  +  ,  . 
dt  ^  ^ 

(7) 

In  the  above  equations,  t  is  time,  (x,y,z)  is  a  right-handed  Cartesian  coordinate 
system  with  x  pointing  toward  shore,  y  alongshore,  and  z  upward.  The  corresponding 
velocity  components  are  (u,v,w),  T  is  temperature,  S  is  salinity,  p  is  density,  and  p  is 
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pressure.  Table  1  provides  a  list  of  other  symbols  found  in  the  model  equations,  as  well  as 
values  of  constants  used  throughout  the  study. 

For  the  finite  differencing,  a  space-staggered  B-scheme  (Arakawa  and  Lamb, 
1977)  is  used  in  the  horizontal.  Batteen  and  Han  (1981)  have  shown  that  this  scheme  is 
appropriate  when  the  grid  spacing  is  approximately  on  the  same  order  as,  or  less  than,  the 
Rossby  radius  of  deformation,  which  meets  the  criteria  of  this  study.  The  horizontal  grid 
spacing  is  14  km  in  the  north-south  direction  and  1 1  km  in  the  east-west  direction,  while 
the  internal  Rossby  radius  of  deformation  is  ~  30  km.  In  the  vertical,  the  10  layers  are 
separated  by  constant  z-levels  of  13, 46,  98,  182,  316,  529,  870,  1416,  2283,  and  3656  m. 
This  spacing  scheme  concentrates  more  on  the  upper,  dynamically  active  part  of  the 
ocean,  above  the  thermocline. 

The  model  domain  (Figure  2)  is  a  rectangular  region  encompassing  the  west  and 
south  coasts  of  Australia,  from  ~  22.5°S  to  40°S  (1792  km  alongshore),  and  from 
~107.5°E  to  122. 5°E  (1408  km  cross-shore).  The  coastal  boundaries  of  the  model  domain 
are  closed,  and  have  both  the  tangential  and  normal  components  of  velocity  set  to  zero. 
Bottom  topography  has  been  omitted  to  focus  on  the  roles  played  by  thermal  and  saline 
forcing.  The  constant  depth  used  in  the  model  is  4500  m.  A  modified  version  of  the 
radiation  boundary  conditions  of  Camerlengo  and  O'Brien  (1980)  is  used  for  the  open 
ocean  domain  boundaries  to  the  north,  south,  east  and  west.  Some  spatial  smoothing  is 
applied  in  the  vicinity  of  the  open  boimdaries. 

The  model  uses  biharmonic  lateral  heat  and  momentum  diffusion  with  the  same 
choice  of  coefficients  (i.e.,  2.0  x  10‘’cm'^s'‘)  as  in  Batteen  et  al.  (1989).  Holland  (1978) 
showed  that  the  highly  scale-selective  biharmonic  diffusion  acts  predominantly  on 
submesoscales,  while  Holland  and  Batteen  (1986)  found  that  baroclinic  mesoscale 
processes  can  be  damped  by  Laplacian  lateral  heat  diffusion.  As  a  result,  the  use  of 
biharmonic  lateral  diffusion  should  allow  mesoscale  eddy  generation  via  barotropic 
(horizontal  shear)  and/or  baroclinic  (vertical  shear)  instability  mechanisms.  As  in  Batteen 

2  j 

et  al.  (1989),  weak  (0.5  cm  s’ )  vertical  eddy  viscosities  and  conductivities  are  used. 
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Bottom  stress  is  parameterized  by  a  simplified  quadratic  drag  law  (Weatherly,  1972),  as 
in  Batteen  et  al.  (1989). 

2.  Method  of  Solution 


Equations  (1)  through  (7)  comprise  a  closed  system  of  seven  scalar  equations  and 
seven  unknowns,  u,  v,  w,  p,  p,  T.  and  5.  The  variables,  u,  v,  T,  and  S  are  prognostic 
variables  whose  time  rates  of  change  are  predicted  from  (1),  (3),  (6)  and  (7),  respectively. 
Although  the  diagnostic  variables  w,  p,  and  p  can  be  determined  from  (3),  (4),  and  (5), 
respectively,  there  are  additional  constraints  imposed  on  p  and  w  by  the  choice  of  the 
rigid  lid  boundary  conditions.  The  vertically  integrated  pressure  can  no  longer  be 
obtained  by  integrating  the  hydrostatic  equation  (4)  for  the  free  surface.  Further,  the 
vertically-integrated  horizontal  velocity  is  constrained  to  be  non-divergent,  i.e.. 


fO  f 


ds=0  , 


(8) 


which  is  obtained  by  integrating  (3)  and  applying  the  vertical  boundary  conditions. 

For  any  quantity  q,  let  its  vertical  average  be  denoted  by  q  and  its  departure 

(vertical  shear)  by  q' .  From  (8)  the  vertical  mean  flow  can  then  be  described  by  a 
streamfunction  y/ ,  such  that: 


1  dy/ 


(9) 


1 

V  = - — 

H  ac 


(10) 


The  streamfunction  y/  is  predicted  from  the  vorticity  equation,  which  is  derived 

from  the  vertical  average  of  (1)  and  (2).  Applying  the  curl  operator  to  the  vertical  average 
(1)  and  (2),  and  using  (9)  and  (10),  the  vorticity  equation  is 
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d  \  ^ ^\j/^  1  r dy/ ^  ^  dy/ cH  ' 
dt  H\dx^ )  j  dx.  Sc  ^  c^ 


where  G  and  F  represent  the  collected  contributions  of  the  nonlinear  and  viscous  terms  in 
equations  (1)  and  (2). 

The  vorticity  equation  (11)  is  solved  by  obtaining  an  updated  value  of  ^  by 

application  of  the  leapfrog  (or  every  1 1  time  steps,  the  Euler-backward)  time-differencing 
scheme.  The  associated  value  of  y/  can  then  be  obtained  from: 

H\dx^)  H\dy^J  dx  dx  dy  dy  ’’ 
which  is  an  elliptic  equation.  A  solution  to  (12)  is  fully  prescribed  by  specifying  the 
values  of  y/  on  the  open  and  closed  boundaries  of  the  model  domain.  Currently,  to  solve 
(12),  the  model  uses  an  elliptic  solver  when  there  are  no  variations  in  coastline  geometry 
and/or  topography,  and  successive  over-relaxation  techniques  when  there  are  variations  in 
coastline  geometry  and/or  topography. 

The  vertical  shear  current  («',  v')  is  predicted  from  (1)  and  (2)  after  subtracting  the 
vertical  mean.  The  results  are: 


du' 

dt 


-14?' 


A  ^ 


•  +  F-F 


dv' 

IF 


r  4  d~V'  - 

+  — +  G-G 


(13) 

(14) 


In  (13)  and  (14),  p',  which  represents  the  departure  of  the  pressure  from  the 
vertical  average,  is ,  using  (4),  expressed  in  terms  of  p  as: 
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(15) 


p'  =  lpgds-^[^{pgds)dz 

where  5  is  a  dummy  variable  representing  the  vertical  coordinate. 

The  method  of  solution  consists  of  predicting  u'  ,v' ,  T,  and  S  from  (1 1), 

(12),  (13),  (14),  (6)  and  (7),  respectively.  The  total  current  is  then  obtained  by  adding  the 
vertical  shear  part  to  the  vertical  average  part,  after  the  latter  is  obtained  from  y/  using 

(9)  and  (10).  The  diagnostics  p,  w,  and  p’  are  obtained  from  (5),  (8),  and  (15), 
respectively. 

3.  Forcing  Conditions  and  Experimental  Design 

Previous  experiments  by  Batteen  and  Rutherford  (1990)  and  Batteen  et  al.  (1992) 
investigated  the  role  of  thermal  forcing  in  generating  the  Leeuwin  Current  and  eddies  off 
Western  Australia  with  a  model  domain  which  was  a  closed  eastern  boundary  with  open 
borders  to  the  north,  south,  and  west.  The  model  was  initialized  with  climatological 
thermal  forcing  from  Levitus  (1982),  and  the  ocean  was  then  allowed  to  geostrophically 
adjust  in  the  absence  of  external  forcing. 

In  a  more  recent  study  by  Batteen  et  al.  (1996),  thermal  forcing  conditions  from 
Levitus  et  al.  (1994)  and  Levitus  and  Boyer  (1994)  were  used.  In  addition,  the  forcing 
was  continuously  applied  at  the  western  boundary,  and  the  model  domain  was  expanded 
to  include  the  Southern  Australian  region,  an  area  rich  with  observations  of  currents  and 
eddies,  but  lacking  in  eddy-resolving  modeling  studies  of  the  region.  Model 
developments  included  the  opening  up  of  the  eastern  boundary  in  the  vicinity  of 
Esperance  in  Southern  Australia,  and  the  incorporation  of  an  irregular  coastline  into  the 
model. 

In  this  study,  to  explore  the  additional  effect  of  salinity  on  the  LCS,  both  annual 
and  seasonal  salinity  forcing  conditions  from  Levitus  et  al.  (1994)  and  Levitus  and  Boyer 
(1994)  have  been  incorporated  into  the  model.  The  annual  (seasonal)  temperature  and 
salinity  forcing  conditions  for  the  upper  five  levels  (top  level),  which  are  assumed  to  be 
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zonally  homogenous,  are  shown  in  Figure  10  (Figure  11).  For  both  annual  and  seasonal 
forcing  conditions,  since  the  lower  five  levels  do  not  exhibit  much  horizontal  variation, 
they  are  assumed  to  be  constant  for  each  level.  The  temperature  values  used  for  levels  6 
to  10  are  9.53°C,  6.03'’C,  3.24°C,  2.19°C,  and  1. 27*^0.  The  salinity  constant  used  for  the 
lower  five  levels  is  34.7. 

The  solid,  dashed,  and  dotted  lines  plotted  in  Figure  1 1  correspond  to  latitudinal 
regions  where  the  three  water  masses  discussed  earlier  are  found.  Regardless  of  the 
region,  the  temperature  conditions  (Figure  11a)  indicate  that  there  is,  as  expected,  a 
temperature  maximum  in  March  (austral  summer)  and  a  temperature  minimum  in 
September  (austral  winter)  throughout  the  whole  region.  ITie  salinity  conditions  (Figure 
1  lb)  show  relatively  strong  aperiodic  (non-seasonal)  variations  throughout  the  year  for  all 
of  the  regions  presented. 

The  design  of  the  model  experiments  (see  Table  2)  is  as  follows.  Experiments  1 
and  2  study  the  model  response  to  annual  forcing.  Experiment  1  uses  the  variable  thermal 
and  salinity  forcing  conditions  shown  in  Figures  10a  and  10b.  Experiment  2  uses  the 
same  variable  thermal  forcing  conditions  but  has  salinity  set  to  a  constant  such  that  S  is 
replaced  by  Sq  in  equation  (5).  Experiments  3  and  4  study  the  model  response  to  seasonal 
forcing.  Experiment  3  uses  the  variable  thermal  and  salinity  forcing  conditions  shown  in 
Figures  1  la  and  1  lb.  Experiment  4  uses  the  same  variable  thermal  forcing  conditions  but 
has  salinity  set  to  a  constant  such  that  S  is  replaced  by  Sq  in  equation  (5). 

B.  RESULTS  FROM  MODEL  SIMULATIONS 

In  each  of  the  experiments,  the  model  integrations  start  from  a  state  of  rest.  Once 
a  day,  the  model  is  updated  at  the  western  boundary  with  annual  or  seasonal  forcing 
conditions.  To  ensure  that  the  time  period  when  the  Leeuwin  Current  is  strongest  (~  April 
to  September),  is  incorporated  in  the  model,  each  model  integration  is  run  for  a  full  year. 

Although  we  are  interested  in  the  difference  between  model  runs  with  variable 
salinity  and  with  constant  salinity,  we  first  show  that  the  model  produces  reasonable 
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results.  The  main  phenomena  in  the  LCS  is  the  density-driven  Leeuwin  Current,  which 
flows  poleward  off  Western  Australia  and  eastward  off  Southern  Australia;  an 
equatorward  undercurrent  found  below  the  Leeuwin  Current  off  Western  Australia;  and 
eddies.  Figures  12  and  13  show  that  the  Leeuwin  Current  with  velocities  of  -  30-100 
cm/s,  an  equatorward  undercurrent  with  core  velocities  of  -  5-15  cm/s,  and  anticyclonic 
eddies  (one  just  south  of  Shark  Bay,  one  between  Dongara  and  Fremantle,  and  one  in  the 
vicinity  of  Cape  Leeuwin)  are  reproduced  in  the  model.  The  results  from  model  runs  are 
qualitatively  similar  and  compare  well  with  available  observations  of  the  LCS  (e.g., 
Boland  etal,  1988;  Cresswell  and  Golding,  1980). 

1.  Experiments  1  and  2  -  Annual  Thermal  and  Salinity  Forcing 

Experiments 

The  annual  forcing  model  integration  with  variable  salinity  (Experiment  1)  and 
that  with  constant  salinity  (Experiment  2)  begin  with  identical  temperature  distributions. 
When  the  two  integrations  are  compared  with  each  other,  significant  differences  are 
found,  as  is  illustrated  here  for  the  time  period  that  the  Leeuwin  Current  is  strongest  (  - 
April  to  September,  corresponding  to  model  days  90  to  273).  Differences  of  up  to  4°C  are 
noted  during  this  time  period  (Figure  14a),  due  to  the  advection  of  warmer  (21  °C  or 
greater)  water  in  Experiment  1  by  the  Leeuwin  Current  between  Cliffy  Head  and  Albany 
(compare  Figure  15a  with  15b).  Farther  east,  between  Albany  and  Esperance,  the 
temperature  difference  of  —  3°C  in  Experiment  1  is  due  to  the  trapping  of  wanner  water 
nearshore  by  the  Leeuwin  Current. 

The  time-averaged  salinity  field  (Figure  16)  shows  that  the  three  major  water 
masses  discussed  earlier  (see  Figure  2)  are  reproduced  well  in  the  LCS  model,  i.e.,  fresh 
(<35.5)  water  is  found  off  the  northern  part  of  Western  Australia,  salty  (>35.5)  water 
offshore  of  Western  Australia  and  Cape  Leeuwin,  and  fresh  (<35.5)  water  to  the  south. 

The  time-averaged  density  fields  (Figures  17a  and  17b),  along  with  the  density 
difference  between  Experiments  1  and  2  (Figure  14b)  show  the  following.  The  density 
differences  north  of  34°S  are  relatively  small.  The  water  for  Experiment  1  off  the 
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northern  part  of  Western  Australia  is  less  dense  than  that  in  Experiment  2,  due  to  the 
combination  of  warm  and  fresh  water  in  the  region,  while  offshore  of  Western  Australia, 
the  water  for  Experiment  1  is  more  dense,  due  to  the  combination  of  salty  and  warm 
water. 

The  largest  density  differences  (Figure  14b)  are  found  between  Cape  Leeuwin  and 
Esperance.  The  water  off  Cape  Leeuwin  in  Experiment  1  is  more  dense  than  that  in 
Experiment  2  due  to  the  saltier  water.  The  water  between  Cliffy  Head  and  Albany  is  less 
dense  in  Experiment  1  than  in  Experiment  2,  due  to  the  much  warmer  (up  to  4°C)  water 
in  the  region.  The  water  east  of  Albany  is  more  dense  in  Experiment  1  than  in 
Experiment  2,  due  to  the  saltier  water  nearshore  which  has  been  advected  into  the  region 
by  the  Leeuwin  Current. 

2.  Experiments  3  and  4  -  Seasonal  Thermal  and  Salinity  Forcing 

Experiments 

Since  the  LCS  has  significant  temporal  variability  in  salinity  (see  Figure  lib). 
Experiments  3  and  4  are  designed  to  study  the  model  response  to  seasonal  forcing,  both 
using  the  same  thermal  forcing  conditions  but  with  variable  (Experiment  3)  and  constant 
(Experiment  4)  salinity.  As  in  the  previous  experiments,  here  we  focus  on  the  difference 
between  the  two  model  runs  to  identify  the  salinity  effect  on  density.  During  the  time 
period  that  the  Leeuwin  Current  is  strongest,  differences  of  up  to  2°C  are  noted  (Figure 
18a)  due  to  the  displacement  of  warm,  anticyclonic  eddies  by  the  Leeuwin  Current  in  the 
Southern  Australian  region. 

North  of  34°S,  the  density  differences  (Figure  18b)  show  that  the  water  for 
Experiment  3,  compared  to  Experiment  4,  is  less  dense  off  Western  Australia  and  more 
dense  offshore  of  Western  Australia.  The  water  off  Cape  Leeuwin  in  Experiment  3  is 
more  dense  than  that  in  Experiment  4  due  to  the  combination  of  salty  and  warm  water  in 
the  region.  The  water  between  Cliffy  Head  and  Albany  is  less  dense,  due  to  the 
displacement  of  warm  (up  to  2°C),  anticyclonic  eddies  by  the  Leeuwin  Current  in  the 
region. 


16 


IV.  SUMMARY 


This  study  used  climatological  temperature  and  salinity  fields  to  calculate  the 
salinity  contribution  to  density  and  dynamic  height  fields  in  the  Leeuwin  Current  System 
(LCS).  While  the  temperature  gradient  is  primarily  linear,  with  warmest  water  to  the 
north,  the  salinity  field  is  spatially  inhomogeneous  with  low  salinity  (tropical)  water  to 
the  north  and  northwest,  high  salinity  (subtropical)  water  to  the  southwest,  and  low 
salinity  (sub-Antarctic)  water  to  the  south.  A  comparison  of  annual  density  fields, 
calculated  with  constant  and  variable  salinity,  showed  that,  off  Western  Australia,  the 
density  field  was  primarily  determined  by  temperature.  South  of  34*^8,  the  effect  of 
salinity  on  density  became  significant.  In  particular,  off  Southern  Australia,  the  density 
was  shown  to  be  influenced  by  both  the  warm  and  salty  (subtropical)  water  mass  and  the 
fresh  and  cold  (sub-Antarctic)  water  mass.  The  seasonal  density  field  patterns  were  found 
to  be  similar  to  the  annual  density  field  patterns. 

The  annual  dynamic  height  fields,  calculated  with  constant  and  variable  salinity, 
showed  similar  flow  patterns  off  Western  Australia  but  different  flow  patterns  between 
34°S  and  36°S.  Throughout  the  year,  the  seasonal  dynamic  height  fields  calculated  with 
constant  and  variable  salinity  showed  similar  flow  patterns  off  Western  Australia  and 
different  flow  patterns  south  of  36°S. 

In  addition  to  the  analysis  of  the  climatological  temperature  and  salinity  fields,  a 
high  resolution,  multi-level,  primitive  equation  model  was  used  to  investigate  the  role  of 
salinity  in  the  formation  of  currents  and  eddies  in  the  LCS.  After  incorporating  salinity 
into  the  model,  two  identical  ocean  models,  one  with  a  climatological  salinity  field  and 
the  other  with  no  horizontal  salinity  gradients,  were  run  and  compared  with  each  other. 

Experiments  1  and  2  used  annual  thermal  and  salinity  conditions  to  initialize  and 
force  the  model.  Despite  the  model  runs  being  initialized  with  the  same  temperature 
distributions,  differences  of  up  to  4°C  were  noted  during  the  time  period  that  the  Leeuwin 
Current  was  strongest  in  the  Southern  Australian  region,  due  to  the  advection  of  warmer 


17 


water  by  the  Leeuwin  Current  between  Cliffy  Head  and  Albany.  Large  density 
differences  were  also  found  between  Cape  Leeuwin  and  Esperance,  where  the  salinity 
effect  on  density  became  significant. 

Experiments  3  and  4  used  seasonal  thermal  and  salinity  conditions  to  initialize 
and  force  the  model.  During  the  time  period  that  the  Leeuwin  Current  was  strongest, 
differences  of  up  to  2°C  were  noted  in  the  Southern  Australian  region,  due  to  the 
displacement  of  warm,  anticyclonic  eddies  by  the  Leeuwin  Current  nearshore.  Large 
density  differences  were  also  found  between  Cape  Leeuwin  and  Albany,  where  the 
salinity  effect  on  density  became  significant. 

The  results  from  both  the  climatological  analysis  and  the  model  experiments 
strongly  suggest  that  there  are  important  quantitative  differences  between  constant  and 
variable  salinity  calculations,  particularly  in  the  Cape  Leeuwin  and  Southern  Australian 
regions.  The  low  and  high  salinity  tags  in  the  LCS  are  found  to  be  particularly  useful  for 
identifying  the  contributions  of  the  different  water  masses  to  the  development  of  the 
Leeuwin  Current  and  eddies. 


Figure  1.  Box  shows  the  domain  of  the  model  for  the  Leeuwin  Current  System 
(LCS)  off  Western  and  Southern  Australia.  Domain  is  bounded  by  22.5°S  to  40°S, 
107.5°E  to  122.5°E.  The  irregular  coastline  (solid  line)  used  in  the  model  includes 
Shark  Bay  (SB),  Dongara  (D),  Fremantle  (F),  Cape  Leeuwin  (CL),  Cliffy  Head 
(CH),  Albany  (A),  and  Esperance  (E). 
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Figure  2.  Schematic  chart  of  the  mean  large-scale  circulation  in  the  eastern  Indian 
Ocean.  The  tropical  (salinity  <  34.5)  and  subtropical  (>  35.8)  water  masses  are 
shaded.  The  edge  of  the  continental  shelf  (200  m  isobath)  is  dashed.  The  Leeuwin 
Current  is  the  southward  flow  of  warm  low-salinity  water  down  the  Western 
Australian  coast,  (from  Pearce  and  Cresswell,  1985) 
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DEPTH  :  AOOm 


depth  :  400m 


Figure  3.  Annual  temperature  (a)  and  salinity  (b)  fields  at  400  m  depth.  Annual 
density  field  with  constant  (c)  and  variable  (d)  salinity  at  400  m  depth.  The  contour 
interval  is  0.2°C  in  (a),  0.05  psu  in  (b)  and  0.05  in  (c)  and  (d). 
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Figure  5.  Annual  dynamic  height  field  with  variable  salinity  at  (a)  10,  (b)  100,  (c) 
200,  and  (d)  400  m  depth.  The  contour  interval  is  2  dyn.  cm. 
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Figure  6.  Seasonal  density  (at)  fields  with  constant  (a)  and  variable  (b)  salinity  at  10 
m  depth  in  February,  May,  August,  and  November.  The  contour  interval  is  0.5 
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Figure  8.  Seasonal  dynamic  height  fields  meridionally  averaged  at  108.5°E  from  (1) 
22.5°S  to  29°S;  (2)  30°S  to  35°S;  and  (3)  36°S  to  40°S.  The  dashed  (solid)  line  refers 
to  dynamic  height  calculated  with  constant  (variable)  salinity . 


Figure  9.  Seasonal  dynamic  height  fields  in  the  nearshore  region  at  (1)  113.5“E, 
which  is  latitudinally  averaged  from  22‘’S  to  29®S;  (2)  115.5‘’E,  which  is  latitudinally 
averaged  from  30“S  to  35‘’S;  (3)  and  36.5"S,  which  is  longitudinally  averaged  from 
116®E  to  125®E.  The  dashed  (solid)  line  refers  to  dynamic  height  calculated  with 
constant  (variable)  salinity. 


30 


Latitude  (S) 


Temperature  (C) 


re  10.  The  annual  temperature  (a)  and  salinity  (b)  i 
pper  five  levels.  Since  the  lower  five  levels  do  not  h 
ition,  they  are  assumed  to  be  constant  for  each  leve. 


depth  :  13m 
T  (DAY)  :  90 


24.0°S 


26.0''S 


28.0''S  H 


30.0''S  H' 


Ld 

O 

D 


32.0°S  H 


34.0'’S 


36.0''S 


38.0'’S 


109.0°E 


111.0°E  113.0°E  'i15.0°E  m7.0'’E  119.0°E 

100.  longitude 


121.0°E 


Figure  12.  Experiment  1  surface  plotted  every 
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Figure  14.  Difference  of  time-averaged  temperature  (a)  and  density  (b)  from  April  to 
September  (days  90-273)  between  annual  forcing  experiments  with  variable 

(Experiment  1)  and  constant  (Experiment  2)  salinity.  The  contour  interval  is  PC  in 
(a)  and  0.1  in  (b). 
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Figure  15.  Time-averaged  temperature  fields  from  April  to  September  (days  90-273) 
for  (a)  Experiment  1  and  (b)  Experiment  2.  The  contour  interval  is  1°C. 
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Figure  17.  Time-averaged  density  (aj  fields  and  velocity  vectors  from  April  to 
September  (days  90-273)  for  (a)  Experiment  1  and  (b)  Experiment  2.  The  contour 
interval  is  0.5. 
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Figure  18.  Differenee  of  time-averaged  temperature  (a)  and  density  (cr,)  (b)  fields 
from  April  to  September  (days  90-273)  between  annual  forcing  experiments  with 
variable  (Experiment  3)  and  constant  (Experiment  4)  salinity.  The  contour  interval 
is  1°C  in  (a)  and  0.1  in  (b). 
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Plate  1.  Climatological  annual  mean  temperature  (a)  and  salinity  (b)  fields  at  10 
100,  and  200  m  depth.  Note  that  the  coastline  differentiates  with  depth.  Annual 
density  fields  with  constant  (c)  variable  (d)  salinity  at  10, 100,  and  200  m  depth. 
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Plate  2.  Climatological  seasonal  temperature  (a)  and  salinity  (b)  fields  at  10 
depth  in  February,  May,  August,  and  November. 


LATITUDE 


Table  1.  Values  of  Constants  Used  in  the  Model 


Constant 

Value 

Definition 

T, 

278.2°K 

v.^onstant  Keterence  remperature 

So 

34.7 

Constant  Reference  Salinity 

Po 

1.0276  gm  cm^ 

Density  of  Sea  Water  At  Tq  and  A'o 

a 

2.4  X  (°K)" 

Ihermal  Expansion  Coefficient 

p 

7.5  X  lO'^  ~ 

Saline  Expansion  Coefficient 

K 

10 

Number  of  Levels  In  Vertical 

Ax 

1.1  X  lO^cm 

Cross- Shore  Cirid  Spacing 

Ay 

1.4  X  10*’cm 

Alongshore  (ind  Spacing 

H 

4.5  X  lO"'  cm 

Total  Ocean  Depth 

At 

800  s 

1  ime  Step 

fo 

0.75  X  10^  s" 

J _ _ _ T  ■■ 

Mean  Conolis  Parameter 

g 

980  cm  s 

Acceleration  of  Gravity 

Am 

2  X  10''  cm‘‘  s'* 

Biharmomc  Momentum  Diffusion  Coefficient 

Ah 

2  X  10"  cm**  s'* 

Biharmonic  Heat  Diffusion  Coefficient 

0.5  cm^  s'* 

Vertical  Eddy  Viscosity 

Kh 

0.5  cm^  s'* 

Vertical  Eddy  Conductivity 
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Table  2.  The  Design  of  The  Experiments 


Experiment  # 

Experiment 

Response  To  Annual  Forcing 

1 

Variable  Thermal  Forcing  with  Constant  Salinity  Conditions 

2 

Variable  Thermal  and  Salinity  Forcing  Conditions 

Response  To  Seasonal  Forcing 

3 

Variable  Thermal  Forcing  With  Constant  Salinity  Conditions 

4 

Variable  Thermal  and  Salinity  Forcing  Conditions 
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